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Abstract 


In recent years, the significance of oxidative stress in the pathophysiology of 
N eurodegenerative/ developmental disorders like Attention Deficit 
Hyperactivity Disorder, Parkinson’s and Alzheimer’s is being studied at an 
accelerating pace. Nrf2 activation via Keap1 inhibition is an established strategy 
for improving the activity of the cellular antioxidant mechanism. In this study, 
pharmacophore modeling was employed to design efficient Keapl inhibitors 
from well-known polypharmacological phytochemicals after extensive structural 
Article Info modifications to improve their pharmacodynamic, pharmacokinetic and drug- 
Se a te likeness qualities (BBB > 0.9, HIA > 0.85). Density functional theory-based 
Volume 3, Issue 2, April 2021 quantum chemical calculations at the B3LYP/ 6-31G (d, p) level of theory were 
Received : 12 September 2020 performed for the geometry optimization of the novel ligands and for computing 
Accepted : 09 January 2021 their electronic properties. Resveratrol-4 was found to be the most desirable 
Published: 05 April 2021 candidate with and AE = 4.24497 eV. HOMO and LUMO distribution of the 
doi: 10,33472/A FJBS.3.2.2021.130-180 Resveratrol-4 was found to be very favorable for keapl1 binding. Molecular 
docking studies and comparative interaction analysis also ranked the Resveratrol- 
4 derivative as the best multi-domain antagonist of the Keap1 protein with a 
binding affinity of -8 kcal/ mole. The following study presents the application of 
Resveratrol-4 a novel, modified, phytochemical derivative, as an efficient 
antagonist of the Keap1 protein for enhancing nrf2 mediated neuroprotection 
from redox insults. 
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1. Introduction 


Attention Deficit H yperactivity Disorder or as described by Dr. Russell Barkley as the“ Diabetes of the human 
brain,” isaseriously underrated neurological disorder, if untreated, has lifelong socioeconomic repercussions 
likelower self-regulation for deferred gratification, ergo, higher selfishness, and lower stability in relationships, 
increased impulsivity, poor executive and cognitivefunction, poor decision-making ability, lower levels of 
general happiness and quality of living (Barkley, 2014). Etiological studies suggest that Attention Deficit 
Hyperactivity Disorder is physiologically aggravated dueto stress (oxidativeand otherwise caused dueto 
lack of physical/ mental wellbeing and maintenance) among several other factors. 
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Oxidativestress is one of thekey interdeoendent factors that cause neurodegeneration. The combinatorial 
effects of nitrosative and oxidativestress caused by the dysregulation of theneuronal antioxidative mechanism 
can cause various biochemical aberrations that eventually lead to neuronal cytotoxicity and death. Under 
excessive oxidative stress from Reactive oxygen species (ROS) and Reactive nitrogen species (RNS), the 
catecholaminergic catabolism proceeds in undesirable secondary oxidative pathways that culminatein the 
deposition of neurotoxic by-products, by-products which accumulate in the neuroplasm and causes its death 
and eventual necrosis. Some of the neurotoxins produced from such aberrant pathways are dopamine-o- 
quinone, dopaminochrome, and 6-hydroxydopamine from dopamine and norepinephrine-o-quinone, 
norepinephrinesemiquinone, and 6-nitronorepinephrinefrom norepinephrine. Further, the neurotoxic by- 
products from the secondary oxidative pathways with epinephrineas thesubstrate resemblenorepinephrine 
by-products (Napolitano et al., 1999a, 1999b and 2011). 


Additionally, other substrates for oxidative stress-induced catabolism are polyunsaturated fatty acids 
(PUFA) and purines. Unchecked redox imbalanceresults in the excessive catabolism of cerebral docosahexaenoic 
acid (Lauritzen et al., 2016; Hashimoto et al., 2017), which terminates in the production of copious amounts of 
the nontoxic expiratory Attention Deficit Hyperactivity Disorder biomarker, ethane (Ross ef a!., 2003; and 
Yakubenko and Byzova, 2017). Likewise, nucleic acid bases (A denine and Guanine) under similar conditions 
produce8-oxo- 7, 8-dihydroguanine (8-oxoG), or 8-oxodeoxyguanosine (8-oxodG), leading to neuronal DNA 
damage (which could be the etiological basis for the genetic component of Attention Deficit H yperactivity 
Disorder) (Chovanovaé al., 2006). The accumulation of these toxins results in neuronal call death and thereby 
contributes to the pathophysiology and progression of Attention Deficit Hyperactivity Disorder. 


The Nuclear factor erythroid 2-related factor 2 (Nrf2) pathway is the parent cascade that regulates the 
proteins and enzymes that control neuronal antioxidation. The plethora of products of the N rf2 associated 
genes that aid in theneuronal redox homeostasis makes theN rf2 genecascadethe master regulator of neuronal 
redox equilibrium. A recent study on Protandim mediated N rf2 activation in erythrocytes reports a 34% increase 
inthe amount of superoxidedismutase produced (N elson & al., 2006). Previous docking studies conducted for 
nrf2 activation via keap1 inhibition using phytochemicals by Li e¢ al. report that antioxidants are the most 
efficient class of phytochemicals with the potential to inhibit keap1 activity (Li & al., 2019). Current research on 
N rf2 activation by Bhakkiyalakshmi et al. has championed pterostilbene, which was extensively studied tn its 
as-is state for its keap1 inhibition activity (Bhakkiyalakshmi et al., 2016). However, native/ unmodified 
oterostilbene has two reactive methoxy groups, which are highly susceptible to degradation and would 
eventually posethe problem for which it was used to treat (Increasein ROSand RNS concentration). Therefore, 
in the present investigation, Resveratrol (aless reactive and naturally occurring derivative of pterostilbene) 
and other antioxidant phytochemicals were subjected to structural modification and derivatization for 
increasing their binding efficiency to Keap1 activesites and for degradativesef-protection from redox insults. 


Structural modification for drug design has been carried out extensively by reourposing existing drugs 
(Uzzaman and Mahmud 2020), phytochemicals (Gordaliza 2007), peptides (Di 2015), sIRNA (Shukla et al. 
2010), among other biomolecules as therapeutics. Pharmacophore modeling of pre-existent bio-chemicals has 
significant advantages over their traditional ab initio counterparts. It requires fewer iterations of computation 
and the sheer abundanceof information available regarding the basal molecule, likeits potency, drug-likeness 
properties, toxicity, pharmacodynamics, and mode of action. Further, quantum chemical studies based on 
density functional theory arean excellent computational chemistry techniquefor studying the environment 
(solvent, cytoplasm, and gas.) specific structureand electrochemical properties of novel drugs. DFT studies 
provide ample insight into the charge localization, molecular electrostatic potential, reactivity, ionization 
potential, and electron affinity of theligands. These quantum chemical descriptors are highly beneficial for 
understanding the ligands’ redox chemistry during biomolecular interactions (ligand-receptor) and are also 
beneficial in comprehending theligand’s behavior in solvated micro/ macro environments. 


Herein, structurally modified, novel, antioxidant, phytochemical derivatives that induce N rf2 activation 
by inhibiting keap1 activity are presented after a comprehensiveinvestigation by quantum chemical techniques 
and molecular modeling studies that are augmented by drug-likeness screening and pharmacokinetic analysis. 


2. Materials and methods 


2.1, Library generation and structural modification 


A ligand library, comprising of four parent phytochemicals known for their therapeutic and, most importantly, 
antioxidant properties, namely Decussatin (Carsten e& al., 2012), Loliolide (Y ang & al., 2011), Resveratrol 
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(stilbene derived polyphenolic phytoalexin found in grapes and raisins) (Gulcin, 2010) and Scopoletin (Shaw 
et al., 2003), were selected for further investigation. The structures, physical and chemical parameters of the 
phytochemicals were obtained from the PubChem database (Kim & al., 2019). Further, ADMET-SAR (Cheng & 
al., 2012) softwaretool and the CB-ligand platform (Liu et al., 2014) wereemployed to predict the blood-brain 
barrier permeability index of these four phytochemicals. The blood-brain barrier permeability index of the 
phytochemicals was predicted to bestrongly positive by both of thesetools. Therefore, thesewere finalized as 
the parent phytochemicals for further ligand library generation via structural modification. Structural 
modifications were performed primarily using ChemDraw Ultra 12.0(2D) (Cousins, 2011). To increaseefficacy, 
antioxidant activity, and binding efficiency to theactivesites of Keap1, the hydroxy! functional groups present 
inthe parent molecules weretargeted for modification via replacement as they can contributeto post catabolic 
oxidative stress. 


Further, thehydroxyl groups wereserially replaced with trifluoromethyl (-CF,) groups (Muller et al ., 2007; 
Hagmann, 2008; and Purser & al., 2008) and amino (-NH.,) groups. It is also noteworthy to know in this context 
that commercially approved selectiveserotonin reuptake inhibitors (antidepressants), Fluvoxamine (Brand 
name: Luvox) (\VWildeet al., 1993), and Fluoxetine (Brand name: Prozac) (A!ltamura e a!., 1994) contains one 
trifluoromethyl moieties each. Further, itis widely known that; the hydroxyl groups contributeto excess redox 
imbalanceand oxidative damagecompared to nitrogen groups (H ybertson & al., 2011). H ence, trifluoromethyl 
groups (widey known as bio-isosteres) wereused in their placeto increase binding efficiency as the fluorine 
groups would increasethe reach of theleads in the keap1 binding pockets (M eanwe!, 2011). Themodified 2D 
structures were converted into a 3-dimensional format using ChensD Pro 12.0, and thefinal structures were 
saved in .sdf and .pdbfileformats. M odifications to the Decussatin molecule was not performed as none of its 
derivatives passed the blood-brain barrier threshold index as set by the two previously mentioned tools. A 
total of eleven derivatives weredesigned with 15ligands in the final ligand library (Figure 1). Theseligands 
weresubjected to further studies. 


2.2. Druglikeness assessment: AD ME 


TheADMET (Absorption, Distribution, Metabolism, Excretion, and Toxicity) assessment of the prepared library 
was conducted using AdmetSA R (Cheng ef al., 2012) tool. A preliminary screening was conducted based on 
Lipinski’s ruleof five(Lipinski, 2000). Further, the blood-brain barrier permeability parameter was prioritized 
as the most stringent rule (BBB >0.80). Additional parameters like Human intestinal absorption, Caco-2 
permeability, P-glycoprotein substrate and inhibitor, Subcellular localization index, Renal organic cation 
transporter, HERG inhibition, AM EStoxicity, and carcinogenicity were computed for a pervasive assessment 
of pharmacokinetic desirability. Ligands with severe outliers weresimply discarded from thelibrary, and the 
ligand derivatives with outliers that could merit minor structural modifications weresubjected to thesame. 


2.3. Quantum chemical calculations 


Quantum chenical calculations of the drug-likeness screened library wereperformed using theGeneral Atomic 
and Molecular Electronic Structure Systen (GAM ESS-US) software(Gordon and Schmidt, 2005). TheA vogadro 
molecule editor and visualizer software (H anwell et al., 2012) was used to generate input simulation files for 
the GAMESS program. M acM oIPIt software (Bodeand Gordon, 1998) was utilized for visualizing the output 
(.log) files. Initial ligand optimization was performed using the A vogadro software. The system’s pH was 
adjusted to that of cerebrospinal fluid, i.e, 7.2 (Albrecht & al., 2020). The Merck molecular force field 94s 
(MM FF94s) was set up for primary optimization of the ligand geometry with water as the solvent system. 
Secondary optimization using GA MESS was performed using the unrestricted density functional theory with 
Becke’s (B) three-parameter hybrid model (Becke, 1988) along with the Lee, Yang, and Parr’s (LY P) non-local 
correlation functional (B3LY P) (Leee al., 1988) along with Pople 6-31G (d, p) asthe basis set for the geometry 
optimization and electronic property calculation of all ligands. The energies of the frontier molecular orbitals, 
highest occupied molecular orbital (HOMO), lowest unoccupied molecular orbital (LUMO), and bandgap 
werealso calculated at theB3LYP/ 6-31G (d, p) levd of theory. Electronic properties likeedectronic nergy, free 
energy, dipole moment, molecular electrostatic potential, and bond angles and distances were calculated. The 
pictorial illustrations of theelectronic spin density distributions of the HOMO and LUMO represented ligand 
structures weregenerated and color-coded to differentiatethe positiveand negative dectrostatic potentials by 
using the M acM oIPIt visualizer. The Parr and Pearson’s interpretation (Parr and Pearson, 1983; Lee é al., 
1988; and Parr and Yang, 1995) along with the Koopman’s theorem (Koopmans, 1934) were utilized for 
calculating thereactivity parameters viz electron affinity (EA), ionization energy (IA), chemical potential (1), 
hardness (77), softness (co), aectronegativity (y), electrophilicity (@), aectron-accepting power (w*) and eectron- 
donor power (a) of theligands. 
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Figure 1: 2D structures of the ligand library: (A) Decussatin, (B-G ): Loliolide (L, L1-6), (H-L): Resveratrol 
(R, R1-4), (M-O): Scopoletin (S, S1-2), 
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2.4. Preparation of ligands 


Ligand preparation was performed using the UCSF Chimera 1.13.1 suite (Pettersen & al., 2004). Theligands 
wereminimized viathestructural editor tool with 100 steepest descent steps and ten conjugate gradient steps 
with astep sizeof 0.02A. Missing hydrogens and Gasteiger partial charges wereadded. AMBER forcefield 
14SB (AM BERFF14SB) was set up during the whole procedure. Thestructures weresaved in pdbat format for 
further studies. 


2.5. R eceptor selection and analysis of activesite 


TheKeap1 protein is thiol-rich (25 cysteine residues), making the protein generally amphoteric and charge 
dynamic. The Broad complex, Tramtrack, and Bric a brac (BTB) (50-184 amino acid residues) domain of the 
K eapl protein is crucial for its homo-dimerization, activity, and structural integrity. TheBTB site contains the 
positively charged, argininetriad residues (R380, R415, and R483, which are crucial for therecognition of the 
Nrf2 protein) and the histidine residues (H 432, H436, H553, H575), which makes the BTB site very basic 
(Dinkova-Kostova & al., 2002). Further, the Kdch domain (327-611 amino acid residues) in Keap1 (onedomain 
per monomer) is essential for capturing and binding N rf2 protans. The DLG and ETGE domains (23-82 amino 
acid residues) in theN eh2 region of the N rf2 protein are fundamental for it to bind simultaneously to two of the 
Keach domain’s Keap1 homodimer (J noff et al., 2014; and Canning et al., 2015). The Kelch domain or Cysteine 
pocket contains C319, C434, C489, C583, C613, and C624, among other residues, making it hydrophobic and 
amphoteric (Dinkova-Kostova & al., 2002). Molecular docking analysis was performed with thesetwo domains 
(BTB and Kelch) as the active sites to investigate the prepared library for dual-site binding efficiency. The 
structureof the Broad complex, Tramtrack, and Bricabrac (BTB) domain (PDB Code: 5DA D) (Resolution: 2.61 
A) and theK ach domain (PDB Code: 4L 7B) (Resolution: 2.41 A ) of the Keap1 protein were obtained from the 
Protein Data Bank (RCBS) (Berman, 2000). Ramachandran plot analysis was performed to assess the 
stereochemical suitability and torsional stability of the amino acids of the selected protein receptors via the 
MolProbity server (Davis eé al., 2007) and the PROCH CECK tool (Laskowski ef al., 1996). 


2.6. Preparation of receptor 


‘Dock Prep’ of the UCSF Chimera 1.13.1 suitewas employed to preparethe selected receptors for molecular 
docking simulation. Initially, the co-crystalized, (1S,2R)2[(1S)1[(1,3dioxo2,3dihydrol1H isoindol2y!)methy! ] 
1,2,3,4tetrahydroisoquinolincarbonyl Jcyclohexanelcarboxylic acid co-crystallized ligand from the kelch 
domain of the receptor (4L.7B) and the (6aS, 7S, 10aS)-8-hydroxy4methoxy2, 7,10atri methyl15,6,6a, 7,10, 10a- 
hexahydrobenzo [h] quinazoline-9-carbonitrile (TX6201) ligand from the BTB domain of the keap1 receptor 
(SDAD) wereremoved. Secondly, the Dunbrack rotamer library was used to add theincompleteside chains of 
thereceptors. N ext, missing hydrogens wereadded for the completed sidechains. Further, Gastelger partial 
charges wereadded, and theA MBER forcefidd 14SB (AM BERFF14SB) was set up for the preparation process. 
The prepared receptors were saved in pdbat format and weresubject to further analysis. 


2.7.M olecular docking simulation 


Molecular docking simulation of the prepared ligands into the prepared receptor was carried out using the 
widely accepted AutoDock Vina (Trott and Olson, 2009) software tool using the PyRx version 0.8 software 
(Dallakyan and Olson, 2015). Rigid docking was performed to investigate theligand-receptor binding efficiency 
at the lowest degree of freedom for higher reliability. Water molecules in the systen wereignored and wereset 
to not interfere with thesimulation in any way. Thegrid box of the BTB domain (PDB Code: 5DAD) wasfixed 
at X: Y: Z::-25.4541 : -2.91812: 10.0365 and likewise, the grid box for the Kelch domain (PDB Code: 4L 7B) was 
fixed at X :Y :Z ::-1.50946: 3.16751: -27.1727 coordinates. Thesizeof thegrid box was heuristically adjusted 
during each simulation to completely incorporate the ligand and the domain that is under purview. 
Exhaustiveness of the search was set to 10 and the simulation of all ligands including the co-crystalized 
ligands were conducted. The binding energies (kcal/ mol) of the docked conformations with the highest binding 
energy coupled with thelowest RM SD (upper and lower bound) wereidentified, and the sameconformations 
were exported to Accelrys discovery studio (Accelrys Software, Inc., San Diego, CA, USA) for interaction 
analysis. The ligand-receptor interaction analysis was performed, and docking attributes liketypeof bonds, 
Interacting atoms-residues, the bond length was recorded, and pictorial illustrations were obtained. 
Additionally, LigPlot (Laskowski and Swindells, 2011) program was also used to assess the ligand-receptor 
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interactions for scaffolding theinteraction analysis results obtained from the discovery studio by acting as an 
additional validator. 


3. Results and discussion 


3.1. Druglikeness screening and pharmacokinetic desirability assessment 


Thecomplete ligand library qualified the Lipinski’s fiverule, which was the primary parameter considered 
for drug-likeness screening. Further, the blood-brain barrier permeability was the second parameter that was 
prioritized for drug-likeness screening. Thetop fiveligands with the highest BBB permeability were, R4>L5 
>L1>L >L2. Likewise, thetop fiveligands with thehighest human intestinal absorption wereR4>L5>L > 
L2>L4. Further, Caco2 cell line permeability prediction resulted in thefollowing descending order of penetration 
in ligands, R4>L5 >S1>R1>S2. It is noteworthy to observe that R4, L5, and Sl arethetop three different 
ligands (with three different parent molecules) with the highest pharmacokinetic desirability. The additional 
pharmacokinetic analysis resulted in desirable scores for these three ligands concerning their toxicity, 
metabolism, and distribution. From the results obtained from the conducted drug-likeness screening and 
pharmacokinetics analysis, it can be envisaged that the derivatives, R4, S1, and L5, would successfully permeate 
the systemic circulatory system and also penetrate past the blood-brain barrier in potent quantities. Upon 
analysis of the bioavailability radar representations of the ligands, it was found that R4 had minor outliers 
pertaining to its lipophilicity (XLOGP3; -0.7 <XLOGP3 <-45), unsaturation (Fraction Csp3; 0.25 <(Fraction 
Csp3) <1) and insolubility (LogS (ESOL); 0<LogS (ESOL) <1) values. Likewise, Slhad oneoutlier in its level 
of unsaturation. Further, the values of L5 was well within the desirable (pink) region of the bioavailability 
radar. As obtained from the Swiss-A DME software, the bioavailability radar illustrations of these ligands are 
represented in (Figure 2). Theresults for the complete pharmacokinetic (absorption, distribution, metabolism, 
elimination) analysis for the entire ligand library is tabulated and presented (Table 1). 
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Figure 2: Bioavailability radar illustrations: (A) Loliolide (L5); (B) Resveratrol derivative-4 (R4); and (C) 
Scopoletin (S1) 
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Table 1, Pharmacokinetic analysis of the designed ligand library 


Blood-brain Human Caco-2 HERG AMES Carcino- 
Ligands barrier intestinal permeability -gene toxicity genicity 
absorption inhibition (three-class) 


[es [0 
a a 
cf onne [a[osmr [nen | mie | asm 
ff ae eee oe | 
ee 
fe fm | esse [emo | ose | onmr | vans | one 


0.5968 0.9885 0.5737 0.9232 0.8653 0.7800 


3.2. Quantum chemical calculations: DFT studies 





Although theA ufbau principleis used to estimate dectrons’ allotment into theorbitals of an atom, DFT studies 
help understand thecharge localization and chargemigration in a multi-atomic molecular system for varied 
ionic states. Quantifiable electrochemical properties of the molecule at different ionic states are crucial for 
understanding its reactivity, behavior, binding affinity (if in acomplex), and most certainly, its redox chemistry 
and antioxidant potential. Geometry optimization of the ligands using DFT studies alleviated the steric 
hindrancein the molecules, and the most stable coordinates of theligands wereobtained. DFT studies conducted 
on theconstructed ligand library revealed the highest ed ectron affinity for Loliolide-5 at 1.986 eV, followed by 
Resveratrol-4 at 1.8775 eV. Scopoletin-5 was placed at ninth place among all ligands. Likewise, the electron- 
donating power (-) of Loliolideand its derivatives werethehighest among all the other ligands (L5>L >L2 
>L4>L1>L3). This was followed by Resveratrol-4 with the second-best electron-donating power (@-) at 
108.74 eV (R4>R >R1>R2>R3). Scopoletin and its derivatives displayed thelowest transition energy (AE= 
LUMO-HOMO), indicating the highest bioactivity among therest (S2 >S >S1). Meanwhile, resveratrol-4 was 
ranked at 11" place with AE =4.24497 eV. Murugavel e al., after experimental and computational studies, 
have reported a AE value of 4.6431 eV along with an electron affinity of 1.4438 eV as an excellent energy 
parameter for thar novel human topoisomerase lla inhibitor (lM urugavel & al., 2019). Hence, it can bespeculated 
that the Scopoletin derivatives (S, S1-2) and R4 have the highest bioactivity and the potential for being an 
efficient inhibitor sincelower transition energies and higher electron affinities (Top fiveligands (including L5, 
R4, S1) with electron affinities >1.68710 eV and top deven ligands with AE <4.24497 eV) was observed for the 
ligands designed inthis study. The LUMO was found to beevenly distributed on both the benzenerings of the 
R4 molecular structure. |n contrast, the HOM O orbitals werefound to bemajorly distributed in and around the 
central C3 =C4 doubly unsaturated linker bond (Figure 3(A )). Likewise, the HOM O orbitals of theS1 derivative 
werefound to bemajorly distributed on theheterocyclic (-O) pyran (C10-C12-C11) ring. Whereas, the LUMO 
orbitals werefound to accumulate onto the homo (C7-C6-C14) hexacyclicring, which isin between thepyran- 
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2-onering and the adjunct amine (-N H 2) and methoxy (-OCH 3) functional moieties and to amedium extent 
over the adjacent heterocyclic pyran ring (Figure 3(B)). The color-coded contour maps of the molecular 
electrostatic potentials representing the electronegative and electropositive centers of these ligands were 
generated and analyzed (Figure4(A, B)). Thecontour maps of theremaining ligands arerepresented in (Figure 
(S2)) (see A ppendix). Finally, theR4 and S1 derivatives of Resveratrol and Scopoletin werefound to bethe best 
candidates with the highest potential for both antioxidant and inhibitory binding functions. Hence, these 
molecules weregiven more focus for further studies. The quantum reactivity parameters of theligands are 
tabulated in (Table 2). Theligands’ electronic properties and energy profile, as computed by GAMESS, are 
tabulated in (Table 3). 


Table 2. The quantum reactivity parameters (eV ) computed for the ligand library based on K oopman’s theorem 
and Parr and Pearson’s interpretation (IE: lonization Energy; EA: Electron Affinity; EN: Electronegativity; 
CP: Chemical Potential; EAP: Electron Accepting Power; EDP: Electron Donating Power) 


eer OMO|«LUMO IE CP Hardness | Softness Electro | EAP | EDP 
(eV) | (eV) re (eV) me (7) (o) philicity | (@+) | (o-) 
be ve (eV) (eV) (a) (eV) | (eV) (eV) 
a= etettsfatat =o tale 
2 — 

(L) -73 | -LO 7B 1B 408 -408 30 
a fefeefor fe fafa | | oe [lo 
oa [alae [fe [ee [se [ae [ew | mo [a [ow 
so [ole [slo [oe |o|a| | o | os [os loe 
fae eee 

Ey 


Ez ae ; 200 Ee 


8 eo 
(R) -533 -327 


om fae of af [aoe m 
ate eee a ee 
ee ee eee 


spe R4 2 jor 
6. 


B oz 
-601 | -10 1® 36 1604 


4 ua 6. 


ee tae algal etal ay setae | 


Table 3. Molecular formula, M olecular weight (g/mol), Total Energy, Electronic Energy in Hartree, and dipole 


moment (D ebye) of the studied polypharmacological phytochemicals and its modified derivatives. 


S: Name M olecular M olecular Total energy Electronic Dipole 
No. formula weight (H artree) energy moment 
(g/mol) (H artree) (Debye) 


304.29 | -1068.875956 | -2858.402387 | 3.094964 
Loliolide 196.24 -654.1009001 -1641.739902 8.796275 
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Table 3 (Cont.) 


M olecular M olecular Total Energy Electronic Dipole 
formula weight (H artree) energy moment 
(g/mol) (H artree) (Debye) 


WH, NO, 195.26 -634.2080522 -1620.996895 8.114084 


7.418502 


C,,H,NO, 195.26 -634.233399 -1619.974298 8.868181 
C,,H,, NO, 193.24 -633.0297468 -1581.563645 4.599633 


C,,H,NO, 191.18 -686.35691 -1540.10211 7.633748 


C,,HNO, 191.18 -666.112597 -1516.636141 9.267832 


Figure 3: The color-coded frontier molecular orbitals depicting the Highest Occupied M olecular 
Orbital (HOMO) and the Lowest Unoccupied Molecular Orbital (LUMO): (A) Resveratrol derivative-R 4; 
(B) Scopoletin derivative-S1 





3.3. Protein receptor analysis by Ramachandran plot 


Ramachandran plot analysis was performed to assess the selected receptor proteins’ structural suitability and 
torsional stability (PDB: 4. 7B / 5DAD). Theallowed values of y (Psi) against @ (Phi) angles for a particular 
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Figure 4: The color-coded contour maps illustrating the molecular electrostatic potential, the red color 
depicts electropositive centers, the blue regions represent the neutral potentials, and the green areas 
represent the electronegative centers: (A) Resveratrol derivative-R 4; (B) Scopoletin derivative-S1 





amino acid in aprotein was assessed. Therearenumerous structures of the K eap1 protein deposited into the 
Protein Data Bank. For thesuitability assessment of proteins for further in-silico studies, the Ramachandran 
plot plays acritical role(H ooft & al.; 1997). For the K dch inhabiting Keap1 protein structure (PDB: 4L 7B), 100% 
(290/ 290) of the amino acid residues were reported in the allowed region. Further, 96.2% (279/ 290) amino 
acids werereported to bein the favored area. Likewise, for the BTB domain harboring K eap1 protein structure 
(PDB: 5DAD), 100% (119/ 119) of the amino acid residues werereported to bein the allowed region. Further, 
95.8% (114/ 119) of theamino acids werein thefavored area. H ence, thesetwo structures were considered the 
most suitable, and further studies were conducted. The general case representations of the Ramachandran 
Plots of the individual proteins are depicted (Figure5 (A, B)). 












































Figure 5: General case representations of the Ramachandran plot analysis representations of the selected 


protein receptor structures depicting the quadrant allocations of the constituent amino acid residues, as 
obtained from the M olProbity server: (A) BTB domain (PDB: 5DAD ); (B) Kelch Domain (PDB: 4L7B). 
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3.4. M olecular docking simulation 


Theentireligand library was subjected to molecular docking and was studied for ther binding efficiency with 
both the BTB domain and Kelch domain of the K eap1 protein. An efficient lead should have excellent binding 
affinities with both the proteins. Such aligand would prevent the keap1 protein’s homo-dimerization and 
would inhibit the binding and eventual ubiquitination and proteolytic degradation of the N rf2 protein. Upon 
analysis, it was found that the Resveratrol derivative R4 ranked in the first place (lowest root mean square 
deviation: upper and lower bound) for binding affinities with both the BTB (-7.3 kcal/ mol) and Kelch (-8 
kcal/ mol) domains. Thetop fiveligands concerning ther binding relationships with the BTB domain wereR4 
>R >R1>R2 >R3. Unfortunately, the Scopoletin derivatives were ranked at the 7” placein S >S2 >Sl, 
respectively. Likewise, thetop fiveligands concerning their binding affinities with the Kelch domain wereR4 
>D >L5>R >R3>S1. Thetherapeutic outcomeby inhibitory prevention of theK eap1-N rf2 complex formation 
is relatively more significant for enhancing the neuronal antioxidant mechanism than the prevention of the 
Keap1homodimerization. Additionally, Pharmacokinetic analysis and quantum chemical studies based on 
density functional theory have regarded theR4, S1, and L5 derivatives as the best ligands with morepromise 
for exerting adesirablelevel of therapeutic outcomein the paradigm of this study. Hence, thebinding/ docking 
of R4, L5, and Slwith the BTB and Kelch domains were studied in depth. Theinteraction analysis obtained 
from LigPlot and Accelrys discovery studio, for someligands, yidded deviant results; this may be becausethe 
algorithms employed by this softwareto predict/ identify interactions aredifferent from oneanother. Howeve,, 
an extensive literature review promotes Accelrys discovery studio to be much more efficient compared to 
LigPlot. Therefore, the interaction analysis obtained from thediscovery studio was given more prominence. 
N evertheless, the pictorial representations obtained from LigPlot for the interaction analysis of the ligands 
that are not under the spotlight of this study with the selected Keap1 and BTB receptors are presented in 
(FigureS3(A-O) and FigureS4(A-O) (SeeA ppendix). 


3.5, Interaction analysis: BTB domain 


The Bis(trifluoromethy! benzene moiety of the R4 derivativewas found to be sequestered into the BTB active 
site, which apparently is sandwiched between thealpha helices (H 131-A 143 and K 150-Y 162) (Figure6 (A, B)). 
Also, TheR4 Bis(trifluoromethyl)benzene moiety displayed a Pi donor- hydrogen bond and a hydrophobic 
bond with the Cysteine 151 residue of the BTB active site. In-depth studies conducted on thesignificance of the 
cysteine residues in the K eap1 protein by Saito et al. report that the Cysteine 151 residue is indispensable for 
the keap1 inhibition mediated ‘turn-on’ mechanism of the cellular antioxidant system as it is one of three 
significant cellular oxidative stress sensing cysteine residues (Saito ef al., 2016). However, Huertae€t al. claim 
that the cysteine 151 binding mediated inhibition itself is not that significant compared to the remaining 
residues that line the BTB binding pocket for keap1 inhibition mediated Nrf2 accumulation (Huerta é& al., 
2016). Holding that thought, R4 was found to interact with Valine 155, Histidine 154 and 129, Arginine 135, 
Lysine 131, Tyrosine 85, and Glutamine 86 (Good agreement with LigPlot results (Figure 6 (C)). Further, the 
fluoride substituted heterocyclic ring adjoining the furan-2-one ring of the L5 derivative was also found to 
interact with Cysteine 151 residue with a hydrogen bond. L5 also interacted with Lysine 131land Histidine 
(154, 129) residues of the BTB domain (Figure 7 (A, B)) (Good agreement with LigPlot results (Figure 7 (C)). 
Finally, SL was found to interact with Arginine 135, Lysine 131, Histidine 129, Histidine 154, and most 
importantly, Cysteine 151 with various types of bonds (Figure8 (A, B)) (Good agreenent with LigPlot results 
(Figure8(C)). Thecomplete interaction profile is tabulated in (Table (4)). The remaining pictorial illustrations 
from PyM ol and Discovery studio for the BT B domain Is represented in (FigureS5) (seeA ppendix). 


Table 4: Post dock interaction profile obtained from Accelrys Discovery Studio for the ligands-BTB domain 


interactions 


Binding affinity Distance A Bond Interaction Interacting 
(Kcal/M ol) category type residues 


Decussatin -6.1 3.53981 Hydrogen bond CHB Lys131 
a 3.67302 Hydrogen bond CHB His154 


i 3.9301 Hydrogen bond PDHB Cys151 





Interaction Interacting 


type residues 





Srinidhi / Afr.).Bio.Sc. 3(2) (2021) 130-180 Page 141 of 180 
aan Binding affinity Bond 
(Kcal/M ol) category 
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Loliolide (L1) 


4.145 H ydrophobic 
A Ays151 

4.72689 H ydrophobic 
A Lys131 


Loliolide (L2) 





A Cys151 
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Table 4 (Cont.) 


Binding affinity 


Distance A Bond Interaction Interacting 





5.06961 H ydrophobic 


H ydrophobic RO Lys131 


3.60548 


4.38336 H ydrophobic 


Val 132 


a 


Electrostatic P-C; PDHB His154 
ee 
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es ee 
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H ydrophobic 
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Loliolide (L4) 
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Table 4 (Cont.) 
Binding affinity Distance A Bond 
(Kcal/M ol) category 
3.87561 H ydrophobic 


5.41155 H ydrophobic 
4.57043 H ydrophobic 


4.46708 Electrostatic 
3.7492 H ydrophobic 


3.87561 H ydrophobic 
5.41155 H ydrophobic 


4.57043 H ydrophobic 
3.18767 Hydrogen bond 


4.46708 Electrostatic 
3.7492 H ydrophobic 


5.07597 H ydrophobic 


3.87561 H ydrophobic 
5.41155 H ydrophobic 


rd 4.46708 Electrostatic 


Interaction Interacting 


type residues 





Lys131 


H ydrophobic 
H ydrophobic 
H ydrophobic 
H ydrophobic 
Hydrogen bond 
Electrostatic 
H ydrophobic 
H ydrophobic 
H ydrophobic 
H ydrophobic 
Hydrogen bond 


4.61343 Electrostatic 
3.8035 H ydrophobic 


4.74585 H ydrophobic 
3.86137 H ydrophobic 


Resveratrol ( 


Resveratrol ( 
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Table 4 (Cont.) 
Binding affinity Distance A Bond Interaction Interacting 
(Kcal/M ol) category residues 
5.31935 H ydrophobic 


4.59489 H ydrophobic 
2.91884 Hydrogen bond 


2.86593 Hydrogen bond oo 





Resveratrol (R3) 





Lys131 


3.81372 H ydrophobic ——— Ais129 


3.7706 H ydrophobic ——- 
5.49754 H ydrophobic 
4.5229 H ydrophobic Cys151 


3.12433 Hydrogen bond Arg135 
3.88118 Hydrogen bond PDHB Cys151 


3.74339 H ydrophobic Lys131 


5.0608 H ydrophobic His129 
3.95133 H ydrophobic His129 


4.07152 Hydrophobic His154 


Lys131 
Arg135 


Scopoletin (S) 


4.74865 H ydrophobic Cys151 


Scopoletin (S1) 





[Wena [aa 


3.77804 H ydrophobic Lys131d 
5. | 5.0017 H ydrophobic His129 


3.98552 H ydrophobic —— His129 
4.02594 Hydrophobic DPA His154 


Scopoletin (S2) 
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Table 4 (Cont.) 


Binding affinity 


Distance A Bond Interaction Interacting 
(Kcal/M ol) category type residues 
3.77424 H ydrophobic A Lys131 


5.14951 H ydrophobic A Arg135 
4.74865 Hydrophobic A Ays151 


Resveratrol (R4) -7.3 3.12891 Hydrogen bond GIn86 

3.25628 Hydrogen bond Lys131 

3.22821 Hydrogen bond Arg135 
Pp 





4.01622 Hydrogen bond; 


3.85983 Hydrogen bond PDHB Ays151 


4.19043 Hydrophobic SPS Hisl29 
5.08778 H ydrophobic P-P-TS His154 


P 
P 
P 
H 
H 
H 
P 
4.06628 H ydrophobic A Lys131 
4.10535 Hydrophobic A Arg135 
A 
A 
P 
Pp 
H 
H 
A 


4.18138 H ydrophobic Ays151 


5.3425 H ydrophobic 


5.14967 H ydrophobic aA Val155 
3.97664 H ydrophobic Lys131 
A 


3.12891 Hydrogen bond; Halogen} CHB; X GIn86 


(F) 
(F) 
3.22821 Hydrogen bond; 
Halogen CHB; X(F) Arg135 
(F) 





Loliolide (L5) 


l 
+ 
iN 


Hydrogen bond; 
Halogen CHB; X Hisl29 


3.1168 Hydrogen bond; a 


S 

3.19794 Hydrogen bond His129 
B 
S 





3.19794 
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Binding affinity Interaction Interacting 

(Kcal/M ol) category type residues 
es 
ee 
re Ee ee 
a 
a 
OO 


3.86608 H ydrophobic His129 


3.89742 H ydrophobic Lys131 
3.72478 H ydrophobic His129 


ote: CHB: Carbon Hydrogen bond; PDHB: Pi-Donor Hydrogen bond; PPS: Pi-Pi Stacked; A: Alkyl; PA: Pi-Alkyl; HHB 
Conventional Hydrogen bond; P-HB; P-SB: Pi-Donor Hydrogen bond; Pi-Sulphur; P-C: Pi-Cation; P-C; P-HB: Pi-Cation 
Pi-Donor Hydrogen bond; P-P-TS: Pi-Pi T-shaped; P-S: Pi-Sigma; CHB;X(F): Conventional Hydrogen bond; Haloge 
(Fluorine); X(F): Halogen (Fluorine). 


Arg135(A) 


RE Lys131(A) 
a: e- cp 


CystS1(A) 


Figure 6: The pictorial illustrations for the Resveratrol-4 derivative in the BTB domain (PDB: 5DAD), 


including the surface view of the top-ranked docked pose (A), interaction analysis from Accelrys Discovery 
Studio (B) and the same from LigPlot (C). 
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Figure 7: The pictorial illustrations for the Loliolide-5 derivative in the BTB domain (PDB: 5DAD), including 
the surface view of the top-ranked docked pose (A), interaction analysis from Accelrys Discovery Studio (B) 


and the same from LigPlot (C) 
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Figure 8: The pictorial illustrations for the Scopoletin-1 derivative in the BTB domain (PDB: 5DAD), 
including the surface view of the top-ranked docked pose (A), interaction analysis from Accelrys Discovery 
Studio (B) and the same from LigPlot (C) 
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Table 5: Postdock interaction profile obtained from Accelrys Discovery Studio for the Ligands-K elch domain interactions 








Binding affinity Distance A Bond Interaction Interacting 
(Kcal/M ol) category residues 
3.24218 Hydrogen bond 


3.987 Hydrophobic 


3.95953 H ydrophobic 
3.97514 H ydrophobic 


5.01244 H ydrophobic 
4.96178 H ydrophobic 


5.09825 H ydrophobic 
4.95345 H ydrophobic 


3.0705 Hydrogen bond 
2.98659 Hydrogen bond 


5.1654 H ydrophobic 
5.38432 H ydrophobic 


4.56214 H ydrophobic 
5.49526 H ydrophobic 


Loliolide (L1) 2.86297 Hydrogen bond Tne Ser363 


Decussatin 





Ile461 


Loliolide (L) 


4.0761 H ydrophobic 


Loliolide (L2) 


Loliolide (L3) 
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Table 5 (Cont.) 


Binding affinity 


Distance A Bond Interaction Interacting 
(Kcal/M ol) category type residues 
2.87467 Hydrogen bond Ser363 


3.25252 Hydrogen bond Asn414 
3.03981 Hydrogen bond Ser602 


5.20126 H ydrophobic Ala556 
3.89607 H ydrophobic Ala556 


4.34819 Hydrophobic pk Ala556 
4.26968 H ydrophobic eI Arg415 





Loliolide (L4) 


5.05972 H ydrophobic Arg415 
3.31077 Hydrogen bond Gly364 


2.15642 Hydrogen bond Ser602 
2.55491 Hydrogen bond Asn414 


3.31077 Hydrogen bond Gly364 
2.75642 Hydrogen bond Ser602 
2.55491 Hydrogen bond Pp HHB Asn414 


3.31077 Hydrogen bond; 
Halogen HHB; X Gly364 


; 
; 
; 
H 
Xx 
: 
C 


Loliolide (L5) 





2.55491 Hydrogen bond Asn414 


2.79219 Hydrogen bond Ser363 


1.98037 Hydrogen bond Asn414 


4.21464 Electrostatic PPC Arg415 





Resveratrol (R) 


| | 
o o 
Ul & 


l 
OY 
iN 


HB 
HB 
HB 
A 
A 
A 
A 
A 
HB 
HB 
HB 
HB 
HB 
HB 
; X(F) 
HB 
HB 
(F) 
HB 
HB 
P-C 
P-C 
P-S 
PA 
Resveratrol (R1) HB 
HB 
HB 
HB 
P-C 
P-C 
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Table 5 (Cont.) 


Resveratrol (R2) 
Resveratrol (R3) 
4.6338 


Resveratrol Ro 2.95147 





Bond 
category 


Distance A Interaction 


Binding affinity Interacting 


(Kcal/M ol) residues 


Ala556 


Hydrogen bond | oHHB Ser363 





3.72588 H ydrophobic 


5.47801 H ydrophobic 
2.15698 


2.06752 Hydrogen bond 


cf 
re | ae 


4.07268 Electrostatic 
5.63874 
4.69447 
2.31515 Hydrogen bond 
2.37749 
4.15495 Electrostatic 
5.75362 
H ydrophobic 
Hydrogen bond 


Hydrogen bond 


Ser508 


Hydrogen bond 


Hydrogen bond 





Hydrogen bond 


Hydrogen bond | cH Ser508 
Hydrogen bond poche Gly603 


H 
H 
H 
H 
H 
H 
H 
C 
C 
C 
C 


3.94395 Electrostatic 


Arg415 
Ala556 


3.69051 H ydrophobic 


3.92247 H ydrophobic 
3.72052 H ydrophobic 
4.61022 H ydrophobic 
5.21239 H ydrophobic 
4.54073 H ydrophobic 


4.10525 H ydrophobic 


type 
P-S 
PA 
HB 
HB 
P-C 
PA 
HB 
HB 
P-C 
PA 
HB 
HB 
HB 
HB 
HB 
HB 
HB 
P-C 
PPS 
A 
A 
PA 
PA 
PA 
PA 
PS 
PS 





Scopoletin (S) 3.16699 Hydrogen bond 
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Table 5 (Cont.) 


ane Binding affinity Bond Interaction Interacting 
(Kcal/M ol) category type residues 
Scopoletin (S1) -6.5 Hydrogen bond HHB 
Hydrogen bond 


HHB Asn382 
2.7133 Hydrogen bond 


a 3.63098 Hydrogen bond PDHB Ser363 
~~ 4.09363 Hydrophobic 


PPS Tyr334 
5.10585 H ydrophobic 


PPS Tyr334 
3.97733 H ydrophobic A Ala556 
4.85466 H ydrophobic 


A Arg415 
P-S Tyr5/72 
PPS Tyr 334 


PPS Tyr334 
4.94473 H ydrophobic P-P-TS Tyr572 
5.1015 Hydrophobic A Arg415 


N ote: CHB: Carbon Hydrogen bond; PDHB: Pi-Donor Hydrogen bond; PPS: Pi-Pi Stacked; A: Alkyl; PA: Pi-Alkyl; HHB: 
Conventional Hydrogen bond; P-HB; P-SB: Pi-Donor Hydrogen bond; Pi-Sulphur; P-C: Pi-Cation; P-C; P-HB: Pi- 
Cation; Pi-Donor Hydrogen bond; P-P-TS: Pi-Pi T-shaped; P-S: Pi-Sigma; CHB;X(F): Conventional Hydrogen bond; 
Halogen (Fluorine); X(F): Halogen (Fluorine). 


Scopoletin (S2) -6.2 3.58029 H ydrophobic 
ee 4.47884 Hydrophobic 





Figure 9: The pictorial illustrations for the Resveratrol-4 derivative in the Kelch domain (PDB: 4L7B), 
including the surface view of the top-ranked docked pose (A), interaction analysis from Accelrys discovery 
studio (B) and the same from LigPlot (C). 
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Figure 10: The pictorial illustrations for the Loliolide-5 derivative in the Kelch domain (PDB: 4L7B), 
including the surface view of the top-ranked docked pose (A), interaction analysis from Accelrys discovery 
studio (B) and the same from LigPlot (C). 
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Figure 11: The pictorial illustrations for the Scopoletin-1 derivative in the Kelch domain (PDB: 4L7B), 
including the surface view of the top-ranked docked pose (A), interaction analysis from Accelrys discovery 
studio (B) and the same from LigPlot (C). 
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3.6. Interaction Analysis: Kelch Domain 


The R4 (8kcal/ mol) Bis(trifluoromethyl)benzene moiety displayed numerous interactions with the Kelch 
domain’s activesiteresidues. Noteworthy among them werethevarious types of bondsformed with Arginine 
415 and Serine 508 residues deep in the kelch pore of the keap1. Jnoff & al., while depositing the protein 
structure for the kelch domain of the keap1 protein (PDB: 4L 7B), have presented a comparative interaction 
analysis study between thedocking results and their actual co-crystallized (Nrf2-Keap1 complex) inhibitor in 
thekelch domain of the keap1 protein. In that study, the prominenceof the binding induced inhibition facilitated 
by their novel inhibitor interaction with the Arginine 415 and Serine 508 residues in the kelch pore of the 
Keapl protein is elucidated (Jnoff e& al. 2014). With the strong interactions displayed by the R4 
bis(trifluoromethyl)benzene moiety with theA rginine 415, Serine 508, and thesurrounding residues that line 
the kelch pore, a scenario can be envisaged where the R4 bis(trifluoromethyl)benzene moiety occludes the 
kelch pore, thereby rendering the K eap1 protein incapable of complex formation with the N rf2 protein. The 
reactiveyet amphoterickdch porewould beasuitable microenvironment for the charge stabilized nucleophilic 
and hydrophobic 6 - (ring) systems along with the significantly e ectronegativetrifluoromethy! groups of the 
R-4 derivative. Meanwhile, thetrifluoromethyl-benzene moiety past the double unsaturated hinge region was 
found to interact with Phenylalanine577 and Tyrosine 572 residues of thekelch domain at the brim of the pore 
tunnel (Figure9 (A, B, C)). 


Further, the L5 derivative displayed interactions with kelch pore residues like Serine 602, Glycine 364, 
Asparagine 414, and 382 (Figure 10(A, B, C)). Theimportance of these residues has not been experimentally 
assessed for their binding mediated keap1 inhibition capability. H ence, no definite conclusions weredrawn 
with regard to their individual/ synergistic importance for inhibiting the Keap1 activity. Furthermore, S1 
displayed interactions with Kelch residues like Asparagine 414, Arginine 415, Serine 363, Asparagine 382, 
Tyrosine 334, and Alanine 556. The hydrophobic interactions displayed by Asparagine414 and Arginine 415 
arenoteworthy (Figure 11(A, B, C)). Theinter-residueinteractions that did not warrant inclusion wereexcluded 
whileobtaining pictorial illustrations of theinteraction analysis. Thecompleteinteraction profile is tabulated 
in (Table (5). Theremaining pictorial illustrations from PyMol and Discovery studio for the Kelch domain is 
reoresented in (FigureS6). 


4. Conclusion 


In thisinvestigation, Resveratrol, Decussatin, Loliolide, and Scopoletin weresubjected to in silico derivatization 
by structural modification via sequential replacement of the parent hydroxyl groups with trifluoromethyl and 
amine groups resulting in aligand library. Theentireligand library was subjected to pharmacokinetic analysis 
(ADME) and wasfurther modified to obtain maximum drug-likeness qualities (Blood-brain-barrier permeability 
(BBB >9.0), Human intestinal absorption among others (HIA >0.85)). Quantum chanical calculations based 
on density functional theory were performed to optimizethegeometries of the constructed library. DFT studies 
also revealed theligand’s electronic properties (bond lengths, bond angles, HOMO, LUMO, bandgap, ionization 
potential, and electron affinity) ina water solvent system set to the pH (7.2) of the human cerebrospinal fluid. 
R4 (Resveratrol), L5 (Loliolide), and S1 (Scopoletin) derivatives wereattributed as superior molecules concerning 
their potential as an antioxidant/ free radical scavenger and also for binding mediated inhibition of the Keap1 
protein. Further, molecular docking studies of theligand library ontheK ach and BTB domains of theK eap1 
protein was performed. The post dock comparative interaction analysis suggested remarkable active site 
residue-ligand bond formations based on which theR4 (Resveratrol) derivativewas ranked first (-8 Kcal/ mol) 
as the best multi-domain inhibitor of the Keap1 protein. Overall, it can be concluded that the performed 
pharmacokinetic analysis, quantum chemical calculations, and molecular docking studies on theconstructed 
ligand library have produced potent and efficient phytochemical inhibitors (R4, L5, and S1) that can 
successfully cross the blood-brain barrier and enhance neuronal cytosolic Nrf2 accumulation via Keap1 
inhibition. However, thevalidity of these claims and findings is deoendent on tn vitro/ in vivo experiments, 
along with the efficiency and reliability of the utilized software. 


Further, fluoridetoxicity isan issue posed by thenovel derivatives presented here (although the rationale 
for including trifluoromethyl moieties are extensively described previously in this work), the potency of these 
candidates seems to behigh. The conducted investigation and the following results providestrong primafacie 
evidencethat small quantities of these novel candidates possess high therapeutic potential. Additionally, a 
hybrid prodrug formulation that would rectify the problem posed by fluoridetoxicity can bea futurescope of 
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work. Finally, the R4moleculecan bechenically synthesized by following the catalysis protocols described by 
Furuya e al. (2011), and R4can be fortified along with resveratrol in yogurt, characteristics of which were 
experimentally studied by Sonarthi ef al. (2018). 
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APPENDIX 


Figure Sl: The color-coded contour maps illustrating the Highest occupied molecular 
orbitals; HOMO) and the lowest unoccupied molecular orbital (LUMO) of the entire ligand 
library. 
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Figure S2: The color-coded contour maps illustrating the variation of molecular electrostatic 
potentials over the regions of the ligands in the constructed library. 
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Figure S3: The two-dimensional illustrations of the interaction analysis of the individual ligands 
with the BTB domain of the Keap1 protein as obtained from LigPlot. 


Met147(A) 


i 
yy 
te 
4 are 
oe 
= oe : 
5 3 
Oey . ? 
Bs. 2s 
Mots ahe 
aa a 
ae His129(A) 
vf 
a 
ms, 
se, 
Se * 
7 ’ 


* 


Arg135(A a (A) 


ot 


Decussatin Loliolide 





Srinidhi / Afr.J.Bio.Sc. 3(2) (2021) 130-180 Page 163 of 180 






D Gly148(A) 


. 
of 
a* 


nek * 


Val132(A) 





Arg135(A) 


Loliolide (L1) Loliolide (L2) 


 Wal132(A) 


’ 





His129(A) **.. 





BZ Val155(A is yeu 
; — ae Val155(A) 


His154(A) auton 
Loliolide (L3) Loliolide (L4) 





His129(A) 






c7 
Rr 
RE CysI51(A) 
ook cae C10 
ian 1 
avn” 
ma 
Metl47(A) 


Loliolide (LS) Resveratrol (R) 


Srinidhi / Afr.J.Bio.Sc. 3(2) (2021) 130-180 Page 164 of 180 


| His1$4(A) 








er 


as® 
es™ 
=~,aat 








Met147(A) 
Resveratrol (R1) Resveratrol (R2) 
K “si 7 ‘Ma 
: ; : Lys131(A) 





cD 


*"His129(A) 


‘ 
o* 


” Cys151(A) 
.¢] 





oJ 
Tee 
? 


; His154{A)- 
Val132(A) 


Resveratrol (R3 
i) Resveratrol (R4) 


Srinidhi / Afr.J.Bio.Sc. 3(2) (2021) 130-180 Page 165 of 180 





t 


.° Met147(A) 


His129(A) 





Scopoletin (S) Scopoletin (S1) 


e 
# 
ca 


CA 


O 








Arg135(A) 
CG 


°s 
yee 
oe 


«wate, 
Lys131(A) 


. 
+ 


His stn A 


Scopoletin (S2) 





*% 
His129(A) 


Srinidhi / Afr.J.Bio.Sc. 3(2) (2021) 130-180 Page 166 of 180 


Figure S4: The two-dimensional illustrations of the interaction analysis of the individual ligands 


with the Kelch domain of the Keap1 protein as obtained from LigPlot. 
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Figure S5: The three-dimensional illustrations of the interaction analysis of the individual 
ligands with the BTB domain of the Keap! protein as obtained from the Discovery studio and 
the surface view of the docked ligand in the BTB active site as obtained from PyMol. 
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Figure S6: The three-dimensional illustrations of the interaction analysis of the individual 
ligands with the Kelch domain of the Keap1 protein as obtained from the Discovery studio and 
the surface view of the docked ligand in the BTB active site as obtained from PyMol. 
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